Gluon polarisation in the proton 



Steven D. Bass 

Institute for Theoretical Physics, Universitat Innsbruck, Technikerstrasse 25, Innsbruck, A 6020 Austria 

Andrew Casey and Anthony W. Thomas 
CSSM, School of Chemistry and Physics, University of Adelaide, Adelaide SA 5005, Australia 

We combine heavy-quark renormalisation group arguments with our understanding of the nucleon's 
wavefunction to deduce a bound on the gluon polarisation Ag in the proton. The bound is consistent 
with the values extracted from spin experiments at COMPASS and RHIC. 

PACS numbers: 13.60.Hb, 13.88.+e 



Polarised deep inelastic scattering experiments have 
revealed a small value for the nucleon's flavour-singlet 
axial-charge, (7^''|pDis ~ 0.3, suggesting that the quarks' 
intrinsic spin contributes little of the proton's spin. The 
challenge to understand the spin structure of the pro- 
ton [1 6] has inspired a vast programme of theoretical 
activity and new experiments. Why is the quark spin 
content .g^'^lpDis so small ? How is the spin i of the pro- 
ton built up from the spin and orbital angular momentum 
of the quarks and gluons inside ? 

A major topic of investigation has been the role of 
polarised glue in the nucleon, both in terms of its con- 
tribution to the nucleon's spin and possible supression of 
the nucleon's singlet axial-charge through the QCD axial 
anomaly. Key experiments to measure gluon polarisa- 
tion are COMPASS at CERN and PHENIX and STAR 
at RHIC. In this paper we investigate gluon polarisa- 
tion via the charm-quark axial-charge, matching the re- 
sults of heavy-quark renormalisation group with what we 
know of the proton's wavefunction. We suggest a bound 
|A(7(m^)| < 0.3, which is consistent with the results of 
the present experiments. 

We start by recalling the gi spin sum-rules, which are 
derived from the dispersion relation for polarised photon- 
nucleon scattering and, for deep inelastic scattering, the 
light-cone operator product expansion. At leading twist 
the first moment of the gi spin structure function mea- 
sures a linear combination of the nucleon's scale-invariant 
axial-charges g'^^\ g'"^^ and (?^''|inv plus a possible sub- 
traction constant /3oo in the dispersion relation [1]: 
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Here cns and cs are the non-singlet and singlet Wilson co- 
efficients. In terms of the fiavour dependent axial-charges 
2Ms^Aq — {p, s|g7^75g|p, s) the isovector, octet and sin- 
glet axial charges are g\ ~ Au — Ad, g\ — Au + Ad — 



2As and g^^^ \inv/ E{as) = g^iQ'^) = Au + Ad + As. 
Here E{as) = exp Jj, ^da^ 7(as)//3(as) is a renormal- 
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isation group factor which corrects for the (two loop) 
non-zero anomalous dimension 7(as) of the singlet axial- 
vector current 0, Jfi.5 = "7^75" + dlfj-lbd + s^fj^j^s , 
which is close to one and which goes to one in the limit 

g2 ^ oo; I3{as) = - ^11 - (a?/27r) + ... is the QCD 

beta function and j{as) — /(as/7r)^ + ... where / is 
the number of active flavours. The singlet axial charge, 
5A^|inv, is independent of the renormalisation scale and 
corresponds to ^^''(Q^) evaluated in the limit oo. 
The flavour non-singlet axial-charges are renormalisation 
group invariants. 

The isovector axial-charge is measured independently 



in neutron /3-decays {g 
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1.270 ± 0.003 



and the 



octet axial charge is commonly taken to be the value ex- 
tracted from hyperon /3-decays assuming a 2-parameter 
SU(3) fit (5^^^ = 0.58±0.03 [al). The uncertainty quoted 



for has been a matter of some debate There 
is considerable evidence that SU(3) symmetry may be 
badly broken and some have suggested that the error on 
g^^^ should be as large as 25% [l3|- Indeed, prompted 
by the work of Myhrer and Thomas [l3|, which showed 
that the effect of the one-gluon-exchange hyperflne inter- 
action and the pion cloud [3] of the nucleon was to 
reduce g^^^ calculated in the cloudy bag model to near 
the experimental value, a re-evaluation of these effects on 
■9a ^ 9a^ 9a^ including kaon loops led to the value 
g''^^ = 0.46±0.05 Here the reduction from the SU(3) 
value came primarily from the pion cloud. 

Assuming no twist-two subtraction constant, polarised 
deep inelastic scattering experiments have been inter- 
preted in terms of a small value for the flavour-singlet 
axial-charge: Ipois.Q^^-cx) = 0.33 ± 0.03(stat.) ± 
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On 



0.05(syst.) if one uses the SU(3) value for g'-^' 

the other hand, using the value g]^' = 0.46 ± 0.05 from 
SU(3) breaking, the corresponding experimental value of 

.9a IpDis would increase to g^^^ |pDis = 0.36 ± 0.03 ± 0.05. 
In the naive parton model g^'lpDis is interpreted as the 
fraction of the proton's spin which is carried by the in- 
trinsic spin of its quark and antiquark constituents. 



2 



Historically, the wish to understand the suppression of 
relative to , led to considerable theoretical effort 
to understand the flavour- singlet axial-charge in QCD. 
QCD theoretical analysis leads to the formula [H [l7l - l20j 
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Here A(7partons is the amount of spin carried by polarised 
gluons in the polarised proton (agAg ^ constant as 
^ 00 



171 [18|) and Aqpartons mcasurcs the spin car- 



ried by quarks and antiquarks carrying "soft" transverse 
momentum kf ^ P^, m? where P is a typical gluon vir- 
tuality and m is the light quark mass. The polarised 
gluon term is associated with events in polarised deep in- 
elastic scattering where the hard photon strikes a quark 
or antiquark generated from photon-gluon fusion with 
kt ~ O, [23|. It corresponds to the QCD axial 
anomaly in the flavour-singlet axial-vector current. Coo 
denotes a potential non-perturbative gluon topological 
contribution [H associated with the possible subtraction 
constant in the dispersion relation for gi and possible spin 
contributions at Bjorken 2: = 0, that is outside the range 
of polarised deep inelastic scattering experiments. The 
measured singlet axial-charge is 5^^|pDis — 3a ^ ~ 

In the parton model A^partons is associated with the 
forward matrix of the partially conserved axial-vector 
current J™^ evaluated in the light-cone gauge — 
and corresponds to the quark spin contribution extracted 
from experiments using the JET or AB factorisation 
schemes [lH. For each flavour q, this term and the pos- 
sible topological term Coo are renormalisation group in- 
variants [ij. All of the renormalisation group scale de- 
pendence induced by the anomalous dimension, 7(05), is 
carried by the polarised glue term: 
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where all quantities in this equation are understood to be 
defined in the /-flavour theory. Flavour non-singlet com- 
binations of the Ag are renormalisation group invariant 
so that each flavour evolves at the same rate, including 
heavy-quark contributions (g = c, 6, t) . The growth in 
the gluon polarisation Ag 1/as at large is compen- 
sated by growth with opposite sign in the gluon orbital 
angular momentum. The quark and gluon total angu- 
lar momenta in the infinite scaling limit are given by [22| 
J,M = H3//(16 + 3/)}and JgM = i{16/(16 + 3/)}. 
There is presently a vigorous programme to disentangle 
the different contributions involving experiments in semi- 
inclusive polarised deep inelastic scattering and polarised 
proton-proton collisions [1, [l^l . 

Heavy-quark axial-charges have been studied in the 
context of elastic neutrino-proton scattering [2J, |25| and 
heavy-quark contributions to gi at values above the 



charm production threshold [26|-|29j. Charm production 
in polarised deep inelastic scattering is an important part 
of the COMPASS spin programme at CERN [s^. 

Following Eq.(2) we can write the charm-quark axial- 
charge contribution as 



Ac(Q2) = ACpa 
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where Acpartons corresponds to the forward matrix ele- 
ment of the plus component of the renormalisation group 
invariant charm-quark axial-current with just mass terms 
in the divergence (minus the QCD axial anomaly), viz. 
(c7^75c)con = (c7^75c) - with A;^ the gluonic Chern- 
Simons current, and we neglect any topological contribu- 
tion ^ . For scales ml, Acpartons corresponds to the 
polarised charm contribution one would find in the JET 
or AB factorisation schemes. 

The heavy-quark contributions to the non-singlet neu- 
tral current axial-charge measured in elastic neutrino- 
proton scattering have been calculated to NLO in 



Ref. [2J]. For charm quarks, the relevant electroweak 



doublet contribution (at LO) is 



(Ac- As) 



(/=4) 
inv 
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- 0{l/ml). 
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For the LO contribution this is made up from 
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+ 0{l/mi), 
~^4'Hml)gX:''+0{l/ml). (6) 



Here Acinv — A.c{Q^) evaluated in the limit 
00, where the charm-quark axial-charge contribution is 
2Ms^Ac = (p, s|c7^75c|p, s). The change in Asjnv be- 
tween the 4 and 3 flavour theories in Eq.(6) comes from 
the different number of flavours in E{as) for the 4 and 3 
flavour theories. 



Eq.(6) contains vital information about < f^Ag > in 

L J 00 

Eq.(4) if we know the RG invariant quantity Acpartons- 
Indeed, if the latter were zero and if we ignore the NLO 
evolution associated with the two-loop anomalous dimen- 
sion ^{as), then Eq. ([5]) would imply (at LO): 



A5(^-^)(Q^) = |"^ 



(™c)„(0)|(/=3) 
~~~'9a linv ' \') 



m 



^ Any topological contribution will be associated with some of 
Acpartons being shifted to Bjorken x = 0. In general, topological 
contributions are suppressed by powers of 1/ for heavy-quark 
matrix elements l3lil . 
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or Ag ~ 0.23 when as{Q'^) ~ 0.3. The following discus- 
sion is aimed at assessing the possible size of Acpartons 
plus the NLO evolution associated with 7(cks), and hence 
the error on this value. 

Canonical (anomaly free) heavy-quark contributions 
to the proton wavefunction are, in general, suppressed 
by powers of so we expect Acpartons ~ 0{l/ml). 

The RG invariant quantity Acpartons takes the same 
value at all momentum scales. We may evaluate it us- 
ing quark-hadron duality in a hadronic basis with me- 
son cloud methods [32|. Experimental studies of the 



strange quark content of the nucleon over the last decade 
have given us considerable confidence that both the ma- 
trix elements of the vector and scalar charm quark cur- 
rents (wh ich are anomaly free) in the proton are quite 
small [33|, [34 1. This gives us confidence in estimating 
the polarised charm contribution through its suppression 
relative to the corresponding polarised strangeness -0.01 
[l5| by the factor ~ (mA — toat -I- rriKY /Ami < 0-1 so 
that |Acpartons| < 0.001. QCD 4-flavour evolution and 
Eq.(6) then enables an estimate of Ag at scales relevant 
to experiments at COMPASS and RHIC. 

In perturbative QCD the LO contribution to heavy- 
quark production through polarised photon-gluon fusion 
yields 



dxg\ 



0, g2 » ml 
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where is the heavy-quark mass {h = c,b,t). The 
anomalous —as/lix term is cancelled against the canon- 
ical term when ^ (the typical gluon virtuality in 
the proton) [l9|. If the gluon polarisation were large so 
that — ^Ag made a large contribution to the suppression 

of g"^ , at this order one would find also a compensating 
large canonical polarised charm contribution in the pro- 
ton. To understand this more deeply, we note that the 
result in Eq.(8) follows from the complete expression 




Here m is the mass of the struck quark and is the 
gluon virtuality. We next focus on charm production. 
The first term in Eq.(9) is the QCD anomaly and the sec- 
ond, mass-dependent, canonical term gives Acp^',l}.°"g for a 
gluon "target" with virtuality P^ . Evaluating Eq.(9) for 
» P^ gives the leading term dx g'^ ~ ^ | ^ ' 
- hence the result in Eq.(8). 

It is interesting to understand Eqs.(8-9) in terms of 
deriving the QCD axial-anomaly via Pauli-Villars regu- 
larisation (instead of the usual dimensional regularisation 



derivation used in [19]). The anomaly corresponds to the 
heavy Pauli-Villars "quark" , which will cancel against 
the heavy charm-quark for a charm-quark mass much big- 
ger than gluon virtualities in the problem (there are no 
other mass terms to set the scale) . When the axial- vector 
amplitude is evaluated at two-loop level there will be 
gluon loop momenta between and the ultraviolet cut- 
off scale generating a small scale dependence so that the 
cancellation between canonical heavy-quark and anoma- 
lous polarised glue terms is not exact in full QCD. This 
scale dependence corresponds to the two-loop anomalous 
dimension 7(as) in E{oLa). The result in Eq.(8) was 
previously discussed in Refs. [10, 27] in the context of 
the phenomenology that would follow if there were large 
gluon polarisation in the proton. Non-perturbative eval- 
uation of Acpartons allows US to constrain the size of Ag 
given what we know about charm and strangeness in the 
nucleon's wavefunction. 

There is a further issue that the derivation of Eq.(6) 
involves matching conditions where the spin contribu- 
tions are continuous between the 3 and 4 flavour theories 
at the threshold scale TOc modulo 0(l/m^) corrections, 
which determine a theoretical error for the method. Us- 
ing QCD evolution with the renormalisation group factor 
i?(as), the results in Eq.(6) are equivalent to the leading 
twist term Ac(m^) vanishing at the threshold scale m^ 
modulo 0(1/to^) corrections, viz. Ac(r7i^) — 0{l/ml) 
[13]. The leading Oil/ml) term is estimated using ef- 
fective field theory in Refs. [13, [13, [13 ■ For polarised 
photon-gluon fusion, this is the — leading term 

in the heavy-quark limit of Eq.(9). The heavy charm- 
quark is integrated out at threshold to give the matrix 
element of a gauge invariant gluon operator with dimen- 
sion 5 and the same quantum numbers as the axial- 

Ac{ml) 



vector current, viz. 



Ac{ml) ^ 0{aAml)AlJml) 



O 



0.017 



27 or 



Tak- 



ing this as an estimate of the theoretical error gives 
Ac{ml) = 0±0.017. 

We next combine this number for Ac(m^) with our es- 
timate of the canonical charm contribution ]Ac, 



partons 



< 



0.001 in quadrature to obtain a bound including theoret- 
ical error on the size of the polarised gluon contribution: 
]~fjA.g(m2)| <0.017or 



]Ag(m2)| <0.3 



(10) 



with Q!s(m^) = 0.4. Values at other values of are 
readily obtained with Eq.(3) or OfgAg ^ constant for 
large values of Q^. 



^ These 0(l/m^) terms associated with the full Ac are manifest 
in polarised photon-gluon fusion through the heavy-quark limit 



(7*9) 



ofEq.(9), Jodxg); 
from the model evaluation of Acparto 
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and arc to be distinguished 
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It is interesting to extend this analysis to full 



flavour QCD. The valuesof A^^f'^^ 



were derived in Ref. [2j] to NLO in the heavy-quark 
expansion. Taking just the leading-order contribution 
plus the heavy- quark power correction according to the 
recipe [13, 21\ described above gives Ac^^"^)(to^) — 
-0.006 ± 0.017 and Adf=^\m'i) = -0.009 ± 0.017 for 
polarised charm. For the bottom and top quarks one ob- 



± 0.001 ± 0.017, A6(^=6)(^2) 



tains A6(^=5)(to^) 

-0.003±0.001±0.017and At(TO2) = 0±2 x 10-^±0.017. 
Here the first error comes from the 0{l/m\) mass cor- 
rection for the heaviest quark of c, 6, t. The second error 
comes from the other heavy-quarks with lesser mass as 
we evaluate these heavy-quark contributions in terms of 
the measured value of the light-quark quantity g''A\[L-v'^^- 
These numbers overlap with a zero value for §J Agr in the 
relevant /—flavour theories. The QCD scale dependence 
of fjAg starts with NLO evolution induced by Kodaira's 
two- loop anomalous dimension 7(0;^). The combination 
^ Ag is scale invariant at LO. This means that if we work 
just to LO and Ag vanishes at one scale, it will vanish at 
all scales (in LO approximation). The LO QCD evolution 
equation for gluon orbital angular momentum in the pro- 
ton [22[ then simplifies so that Lg(cxD) = i{16/(16+3/)}. 
In practice, the two-loop anomalous dimension generates 
slow evolution of fjAg. Dividing the finite value of this 
combination at large scales by the small value of as gives 
a finite value for the gluon polarisation Ag, which can 
readily be the same order of magnitude as the gluon total 
angular momentum (or larger with cancellation against a 
correspondingly larger gluon orbital angular momentum 
contribution). 

It is interesting that the value of Ag deduced from 
present experiments COMPASS at CERN and PHENIX 
and STAR at RHIC typically give \Ag\ < 0.4 with 
as ~ 0.3 corresponding to | — 3^A(ji| < 0.06 This ex- 
perimental value is extracted from direct measurements 
of gluon polarisation at COMPASS in the region around 

^gluon 

~ 0.1, NLO QCD motivated fits to inclusive 51 
data taken in the region x > 0.006, and RHIC Spin data 
in the region 0.02 < Xgiuon < 0.4. The theoretical bound, 
Eq.(lO), is consistent with this experimental result. 
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